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FOREWORD 


This  study  was  undertaken  as  part  of  'an  investigationAof  navigation 
within  the  solar  system  by  optical  means.-  The  objective  of  the  inves¬ 
tigation  is  to  evaluate' £he  suitability  of  various  physical  phenomena 
as  sources  of  navigational  information^shd  to-  estimate1  the  accuracy^of" 
navigational  information  obtained  by  various  techniques.  The  work  was 
supported  under  Contract  No.  AF  33(61S)-7413  by  the  Navigation  and 
Guidance  Laboratory,  Aeronautical  Systems  Division,  Air  Research  and 
Development  Command,  United  States  Air  Force.  Hr.  William  Harmon  was 
project  monitor  for  the  Navigation  and  Guidance  Laboratory.  The 
Contractor's  report  number  is  GCA  62-24-A.  This  is  the  final  report 
under  AF  33(616)-7413 . 
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ABSTRACT 


(1)  Availabl^Wnformation  concerning  the  noon,  which  is  likely 
to  be  of  interest  to  cosigners  of  navigational  aids  for  lunar  vehicles, 
is  given.  This^include^^ata  on  phases,  spectral  distribution,  size 
and  figurd^ and  the  intensity  of  Lyman  alpha  reflected  from  the  moon . 


(2)  e  Preliminary  design  and  operation  of  an  automatic  sextant  for 
use  in  interplanetary  vehicles  is  discussed,  v The  sextant  is  capable  a, 

of  measuring  angles,  usually  not  exceeding  5°;%itK'  a  precision  of  0.2'??'“*'’^'*- 

4  /i 

(3)  ^Various  aspects  of  satellite  navigation  by  terrestrial 
occultations  of  stars  are  discussed.  A  condition  for  occultation  in 
terms  of  stellar  coordinates  and  the  orbital  elements  of  the  satellite 
is  derived,  also  an  occultation  equation  relating  the  time  of 
occultation  to  the  orbital  elements  is  derived.  The  expected  frequency 
of  occultations  is  discussed.  A  method  of  **virtual  occultatiori"using 
extinction  of  ^tarlight  in  the  atmosphere  is  discussed.  An  error  of 
approximately  t  km  in  the  satellite  position  measurement  is  anticipated 
due  to  unpredictable  variations  of  atmospheric  density.  Finally  the 
interference  due  to  the  brightness  of  earth's  atmosphere  is  discussed. 
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SECTION  I 

INTRODUCTION 


The  purpose  of  this  investigation  is  to  study  various  physical 
phenomena  which  may  be  suitable  sources  of  navigational  information 
for  space  vehicles  within  the  solar  system,  including  the  earth  sat¬ 
ellites.  The  scope  is  limited  to  phenomena  which  can  be  studied  at 
optical  wavelength  by  means  of  instruments  which  are  totally  contained 
in  the  space  vehicle.  During  the  year  1960-1961,  eleven  technical 
reports  were  prepared  and  published  and  the  results  were  summarized 
in  the  Final  Technical  Documentary  Report  entitled  "Investigation  of 
Physical  Phenomena  for  Space  Navigation"  by  Stubbs  (for  a  complete 
list  see  Bibliography). 

This  investigation  was  continued  during  1961-1962,  and  the  topics 
studied  fall  into  three  categories:  (1)  lunar  photometry  and  Lyman- 
alpha  intensity  near  the  lunar  limb,  (2)  theory  and  design  of  an 
automatic  sextant  for  space  navigation  and  (3)  satellite  navigation  by 
terrestrial  occultations  of  stars.  Detailed  results  of  these  studies 
have  already  been  published  in  the  form  of  six  technical  reports.  The 
purpose  of  this  final  report  is  to  summarize  these  results.  For  details 
as  well  as  for  bibliographical  material  the  reader  is  referred  to  the 
original  technical  reports. 


Manuscript  released  by  the  author  December,  1962  for  publicition  as 
an  ASD  Technical  Documentary  Report. 
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SECTION  2 

LUNAR  PHOTOMETRY 


Levy  (GCA  Tech.  Rep.  62- 10- A)  has  summarized  the  available  Infor¬ 
mation  on  the  phases,  spectral  distribution  and  size  and  figure  of  the 
moon.  Figure  2.1  and  2.2  give  the  total  brightness  imd  the  limb 
brightness  of  the  moon  at  various  phases.  The  changes  in  brightness 
with  the  phase  can  be  seen  to  be  very  rapid.  At  quarter  phase  when 
the  moon  appears  practically  as  a  semi-circle  its  brightness  is  roughly 
1/9  times  the  brightness  of  the  full  moon.  The  brightness  of  any  fea= 
ture  on  the  moon  is  proportional  to 


cos  1  I  (limb) 
cos  i  +  X  cos  e 


(2-1) 


whore  i  and  e  are  the  angles  of  incidence  and  reflection  respectively, 

X  *»  0.225  (1  +  tan2  ~)  and  I  (limb)  is  the  intensity  at  the  limb. 

Even  the  hemisphere  of  the  moon  which  is  in  the  sun's  shadow  is 
faintly  Illuminated  by  the  sunlight  reflected  by  the  earth.  This  is 
particularly  noticeable  near  the  new  moon,  and  it  is  possible  to 
measure  the  lunar  limb  by  this  so-called  "ashen  light".  Optimally 
the  surface  brightness  of  the  earthlit  moon  is  5,100  times  less  than 
that  of  the  Bunlit  full  moon. 

Since  the  terminator  of  the  moon  is  usually  Jagged,  all  measure¬ 
ments  should  preferably  be  made  at  the  limb  and  not  at  the  terminator. 

The  reflectivity  of  the  moon  has  bem  measured  at  several  opti¬ 
cal  wavelengths  and  is  found  to  increase  gradually  from  the  blue  to  the 
red  region  of  the  spectrum.  As  a  result  the  relative  spectral  intensity 
distribution  of  the  moonlight  (Figure  2.3)  is  not  exactly  the  same  as 
for  the  sunlight,  but  resembles  the  latter  fairly  well. 

The  figure  of  the  moon  in  the  past  was  thought  to  be  roughly  a 
tviaxial  ellipsoid,  but  recent  measurements  due  to  Schrutka-Rechen- 
stansn  and  Hopkins  indicate  that  within  the  mean  error  of  the  altitude 
measurements  used,  -1.2  km,  the  best  figure  for  the  moon  is  that  of  a 
sphere  of  radius 


u  =  1738.1  t  0.2  km  . 
o 


(2.2) 
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Superimposed  on  this  sphere  are  local  elevations  and  depressions 
of  up  to  ±5  lea,  with  ±3  !a  quite  coccon.  The  lunar  naria  are  an 
average  of  2  kn  lover  than  the  mean  level,  continents  2  ka  higher - 
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SECTION  3 


LYHAN-ALPHA  AND  THE  LUNAR  LIMB 


Levy  (GCA  Tech.  Rep.  62-13-A)  has  made  rough  calculations  of  the 
Intensity  of  the  solar  Lyman-alpha  radiation  reflected  by  the  moon. 

He  finds  an  upper  limit  for  the  lunar  limb  brightness  in  Lyman-alpha 
-4  -2  -l  -1 

to  be  4  x  10  ergs  cm  sec  ster  .  The  Lyman-alpha  glow  of  the 
night  sky,  according. to  measurements  of  Rupperian  et  al,  is 
-3  -2-1  -1 

2.7  x  10  ergs  cm  sec  ster  in  the  antisolar  direction  and 
slightly  higher  in  other  directions.  Ir  the  Lyman-alpha  the  moon 
will,  therefore,  appear  silhouetted  against  a  brighter  night  sky 
background . 
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SECTION  4 

AN  AUTOMATIC  SEXTANT  FOR  SPACE  NAVIGATION 


The  theory,  preliminary  design  and  operation  of  an  automatic 
sextant  has  been  discussed  by  Stubb3  (GCA  Tech.  Rep.  62- 17 -A).  The 
sextant  is  expected  to  measure  angles,  usually  r.oc  exceeding  5  degrees, 
with  an  accuracy  of  0.2"  (1  microradian).  This  order  of  accuracy  is 
achieved  by  using  a  prism  pair,  capable  of  rotation  around  a  common 
axis.  Large  relative  rotations  of  the  prism,  which  can  be  measured 
with  good  precision  by  the  usual  mechanical  methods,  produce  very 
small  deflections  of  the  light  beam.  This  is  illustrated  in  Figure 
4.1  for  a  prism  of  1.5°  deflection  angle.  In  this  manner  up  to  60 
times  improvement  in  accuracy  can  be  achieved  over  that  of  the  usual 
mechanical  methods.  Figures  4.2  and  4.3  illustrate  the  arrangement  of 
the  optical  components  of  the  sextant,  which  consists  of  two  tele¬ 
scopes  mounted  parallel  to  each  other  on  the  same  rigid  platform. 
Relative  rotations  of  the  two  prism  pairs  bring  the  images  of  the 
star  and  the  planet  at  the  center  of  each  of  the  trihedral  mirrors. 

The  angle  between  the  star  and  the  planet  is  calculated  from  the 
relative  rotations  of  the  prisms. 

The  design  requirements  for  an  overall  accuracy  of  r  microradian 
in  the  measurement  are  summarized  below: 


Deflection  of  each  element  of  prism  pair 
Prism  rotation  controlled  to 
Telescope  axis  controlled  to 
Maximum  angle  between  light  sources 
Minimum  angle  between  light  sources 
Parallelism  of  the  two  telescopes 


3° 

8" 

20" 

5° 

30' 

0.2"  (1  microradian) 


The  objective  of  the  program  is  the  accurate  measurement  of  four 
celestial  angles;  th*  angles  between  a  first  planet  and  each  of  two 
nearby  stars  and  a  second  planet  and  each  of  two  nearby  stars.  Know¬ 
ing  these  four  angles  ar.d  the  angle  between  the  two  stars  it  is 
possible  to  derive  navigational  information.  The  identification  of 
individual  stars  is  not  requ'-ed.  All  that  is  necessary  la  to  know 
the  angle  between  the  two  sts;i-9.  This  information  concerning  nearly 
5,300  stars  will  be  carried  r„  the  memory  of  the  computer.  It  turns 
out  that  the  angles  between  all  possible  combinations  of  stars  are 
not  needed.  It  is  sufficient  to  store  in  the  computer  memory  the 
angles  between  three  times  as  many  pairs  as  the  number  of  stars 
chosen,  namely  about  16,000  angles  This  requirement  can  be  met  easily. 


8 


9 


Photoelectric  Detector 


NOTE!  The  trihedral  mirrors  and 
re  forcing  mirrors  have 
three- fo'd  'symmetry 


Figure  4.2.  Arrangement  of  optical  components. 
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Refocusing 

Mirror 


Figure  4.3,  Arrangement  of  optical  components. 
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SECTION  5 

SATELLITE  NAVIGATION  BY  TERRESTRIAL  OCCULTATION  OF  STARS 


Let  denote  the  impact  parameter  of  a  ray  in  a  parallel  beam 

of  starlight,  which  is  received  by  the  satellite  at  an  instant  of 
time  t.  For  an  undeflected  ray  (no  refraction)  r^  equals  the  distance 

of  closest  approach  of  the  ray  measured  from  the  earth's  center. 

The  ingress  or  egress  in  a  terrestrial  occuitation  of  a  star 
is  characterized  by  the  condition  that  the  impact  parameter  r^  equals 

the  radius  of  the  earth  aQ.  The  signal  for  the  occurrence  of  ingress 

(or  egress)  is  provided  by  the  disappearance  (or  re-appearar.oe)  of 
the  star  it  the  edge  of  the  earth.  This  method,  which  works  well 
for  lunar  occupations,  however,  fails  for  terrestrial  occupations, 
primarily  on  account  of  the  fact  that  'he  star  can  just  as  easily 
disappear  behin*  an  absorbing  cloud,  whose  unknown  height  above  the 
earth's  surface  may  be  as  much  as  18  km,  as  behind  the  true  edge  of 
the  earth. 

This  difficulty  can  be  surmounted  if  we  consider  a  ray  with  a 
known  impact  parameter,  say  (aQ  +  20)  km,  thus  avoiding  the  unpredict¬ 
able  clouds.  For  a  spherically  symmetric  atmosphere  and  satellites 
above  100  km,  it  turns  out  that  deflection  and  extinction' of  the  light 
beam,  produced  in  '■he  atmosphere,  are,  for  all  practical  purposes, 
functions  of  the  Impact  parameter  only.  The  impact  parameters  can 
therefore  be  calculated  from  the  change  in  direction  or  intensity  of 
light.  We  will  call  the  latter  method,  the  method  of  "virtual  occul- 
tations". 

The  following  equation  (GCA  Tech.  Rep.  62-18-A)  holds  t-etween 
the  satellite  orbital  elements,  the  equatorial  coordinates,  O  and  5, 
of  the  star  and  th,  time  of  observation  t. 


1 


fv . (1  +  e  cos  v) 


a(l  -  e  ) 


1 

2 


=  cos  i  cos  b  sin  (a  -  >T)  sin  (v  +  u>) 

+  sin  i  sin  6  sin  (v  +  co) 

+  cos  5  cos  ( a  -  s<)  cos  (v  +  <»>)  (5.1) 
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where  a,  e  and  i  are  respectively  the  semi  major-axis,  eccentricity 
and  inclination  to  the  equatorial  plane  of  the  satellite  orbit,  v, 
called  the  true  anomaly,  is  a  known  function  of  a,  e,  T  (time  of  per¬ 
igee  passage)  and  the  time  of  observation  t.  r^  is  the  impact  para¬ 
meter  of  the  observed  ray.  For  a  true  occultation  r^  equals  the 
earth’s  radius  aQ,  and  for  a  "virtual  occultation"  it  is  derived  from 
changes  in  direction,  intensity,  etc.,  of  the  light  ray. 

Six  observations  of  occultation  time  for  different  stars  of 
known  coordinates  are  required  to  obtain  the  six  orbital  elements,  but 
due  to  numerous  uncertainties  a  much  greater  number  is  recommended. 

In  most  cases  an  approximate  orbit  of  the  satellite  is  known; 
sometimes  it  is  the  pre-determined  orbit  before  the  launch.  The  con¬ 
dition  that  an  occultation  of  a  star,  whose  coordinates  are  a  and  6, 
will  occur  sometime  during  an  orbital  revolution  of  the  satellite  is 
given  by 


(1  -  e  j  (cos  i  sin  5  -  sin  (a  -  ft)  sin  i  cos  5) 


+ 


e  cos  tan 


-1 


sin  1  sin  6  +  sin  .  .  -  ft)  cos  1  cos  & 
cos  (Ct  -  ft)  cos  & 


<  1 


(5.2) 


(for  derivation  see  GCA  Tech.  Rep.  62-18-A).  This  equation  makes  it 
possible  to  calculate  the  stars  which  will  be  occulted  and  the  number 
of  such  stars. 


In  Figure  5.1  the  number  of  stars  required  (and  the  corresponding 
limiting  stellar  magnitudes)  are  plotted  against  the  radius  of  the 
satellite  orbit  for  25  occultations  per  period  (N^j)  and  per  hour 

(N  )  for  a  circular  orbit  and  a  uniform  distribution  of  stars  over 
2o 

the  sky.  This  figure  does  not  take  two  effects  into  account.  First 
the  occultation  observations  must  be  confined  to  the  night  side  of 
the  earth,  and  second,  the  starlight,  at  the  time  of  observation,  is 
reduced  in  intensity,  by  say  50%,  due  to  atmospheric  extinction.  The 
net  effect  is  that  the  scale  giving  the  number  of  stars  in  Figure  5.1 
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should  be  multiplied  by  two  and  1.4  should  be  added  to  the  magnitude 
scale.  These  changes  are  incorporated  in  Table  5.1,  which  gives  the 
number  of  stars  and  tile  corresponding  limiting  magnitudes  for'.circu- 
lar  as  well  as  elliptical  orbits  for  two  values  of  perigee  distances. 

The  interference  due  to  brightness  of  the  earth's  atmosphere  is 
discussed  in  GCA  Tech.  Rep.  62-22- A  from  which  Table  5.2  is  reporduced. 
The  interference  due  to  day  airglow,  noctilucent  clouds  and  Rayleigh 
scattering  of  sunlight  can  be  ignored  if  occultation  observations  are 
confined  to  the  night  side  of  the  earth.  The  interference  due  to 
aurorae  can  be  large  at  times,  but  is  variable.  The  interference  due 
to  night  airglow  and  scattered  moonlight,  particularly  near  full  moon, 
e  both  quite  large  and  should  be  further  investigated.  These  sources 
n£  atmospheric  radiance  can  be  partially  eliminated  by  using  techniques 
in  which  the  atmospheric  radiance  is  subtracted  from  the  total  radiance 
due  to  the  star  and  the  atmosphere.  Nevertheless,  night  airglow  and 
scattered  moonlight  may  be  quite  troublesome  in  cccultation  observations. 

The  method  of  measuring  extinction  of  starlight  to  derive  infor¬ 
mation  regarding  the  impact  parameter  of  the  light  ray  is  preferred 
over  the  method  in  which  the  deflection  of  the  ray,  due  to  refraction, 
is  measured.  This  method  is  discussed  in  GCA  Tech.  Rep.  62-21-A.  The 
wavelength  region  0,3  to  0.55  p  is  chosen  for  extinction  measurements, 
since  the  atmoephere  has  negligible  true  absorption  in  this  region, 
and  extinction  can  be  predicted  rather  well.  Figure  5.2  shows  percen¬ 
tage  extinction  of  the  atmosphere  as  a  function  of  impact  parameter 
for  two  wavelengths  (0.375  and  0.520  p).  Seasonal  and  latitudinal 
variations  of  the  atmospheric  density  have  not  been  considered  in 
these  calculations.  This  and  several  other  improvements,  lasted  in 
the  above-mentioned  technical  report,  are  required.  After  various 
suggested  improvements  have  been  incorporated  in  the  calculations, 
the  major  source  of  error  is  expected  to  be  the  deviance  of  the  actual 
atmospheric  densities  from  the  standard  model  used.  The  error  in 
position  measurement  of  the  satellite  due  to  these  deviations  is  of 
the  order  of  1/2  km. 
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TABLE  5.1 


NUMBER  OF  STARS  REQUIRED  AND  THE  CORRESPONDING  INSTRUMENTAL  LIMITING' 
MAGNITUDE  FOR  25  OCCULTATIONS  PER  PERIOD  (N25>  AND  PER  HOUR  (N*5) 'FOR- 
A  FEW  SELECTED  ORBITS  (INCLUDING  ATMOSPHERIC  EFFECTS) 

Perigee  Height  hp  =  200  km 
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ma.1  or 
axis 
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of 
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Limiting 
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graphic) 
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graphic) 
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76 
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66 
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98 
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Stellar  Irradiance 


Figure  5.2.  Extinction  cue  to  Rayleigh  scattering  for 

various  impact  parameters  and  two  wavelengths. 
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